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ABSTRACT 


A  shear  deformable  finite  element  is  developed  based  on  generalized 
(i.e. ,  account  for  shear  deformations)  first-approximation  shell  theories. 
Numerical  results  are  presented  for  bending  of  layered,  anisotropic, 
composite  shells.  Various  shell  theories  (e.g.,  Sander's,  Love's, 
Donnell's,  etc.)  are  included  as  special  cases  of  the  present  element, 
and  can  be  obtained  by  giving  numerical  values  to  appropriate  tracers. 

The  present  finite-element  solutions  are  compared  with  the  exact  solu¬ 
tions  for  certain  shell  problems,  and  other  solutions  available  in  the 
literature.  The  agreement  is  found  to  be  very  good.  Numerical  results 
showing  the  effect  of  orientation  of  layers,  boundary  conditions,  and 
material  orthotropy  on  deflections  are  presented. 

INTRODUCTION 

Composite  materials  and  reinforced  plastics  are  increasingly  used 
in  automobiles,  space  vehicles,  and  pressure  vessels.  With  the  increased 
use  of  fiber-reinforced  composites  as  structural  elements,  studies 
involving  the  thermomechanical  behavior  of  composite-material  shells  are 
receiving  considerable  attention. 

The  first  analysis  that  incorporated  the  bending-stretching  coupling 
(due  to  unsymmetric  lamination  in  composites)  is  due  to  Ambartsumyan 
[1,2].  In  his  analyses  Ambartsumyan  assumed  that  the  individual  ortho¬ 
tropic  layers  were  oriented  such  that  the  principal  axes  of  material 
symmetry  coincided  with  the  principal  coordinates  of  the  shell  reference 
surface.  Thus,  Ambartsumyan's  v*ork  dealt  with  what  is  now  known  as 
laminated  orthotropic  shells  rather  than  laminated  anisotropic  shells; 
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in  laminated  anisotropic  shells  the  individual  layers  are,  in  gene-al , 
anisotropic  and  the  principal  axes  of  material  symmetry  of  the  indi¬ 
vidual  layers  do  not  coincide  with  the  principal  coordinates  of  the 
shell. 

In  1962  Dong,  Pister  and  Taylor  [3]  formulated  a  theory  of  thin 
shells  laminated  of  anisotropic  material.  The  theory  is  an  extension  of 
the  theory  developed  by  Stavsky  [4]  for  laminated  anisotropic  plates  to 
Donnell's  shallow  shell  theory.  Cheng  and  Ho  [5]  presented  an  analysis 
of  laminated  anisotropic  cylindrical  shells  using  Flugge's  shell  theory. 
A  first  approximation  theory  for  the  unsymmetric  deformation  of  non- 
homogeneous,  anisotropic,  elastic  cylindrical  shells  was  derived  by 
Widera  and  Chung  [6]  by  means  of  the  asymptotic  integration  of  the 
elasticity  equations.  For  a  homogeneous,  isotropic  material,  the 
theory  reduces  to  the  Donnell  equations.  An  exposition  of  various  shell 
theories  can  be  found  in  the  article  by  Bert  [7]. 

All  of  the  theories  discussed  above  are  based  on  Kirchoff-Love's 
hypotheses  in  which  the  transverse  shear  deformation  is  neglected. 

Recent  studies  in  layered  anisotropic  plates  show  (see  the  survey  paper 
by  the  author  [8])  that  transverse  shear  deformation  effects  are  more 
pronounced  in  composite  plates  and  shells  than  in  isotropic  plates  and 
shells.  The  effect  of  transverse  shear  deformation  and  transverse 
isotropy,  as  well  as  thermal  expansion  through  the  shell  thickness  were 
considered  by  Zukas  and  Vinson  [9]  and  Dong  and  Tso  [10].  The  theory  in 
[10]  is  only  applicable  to  layered,  orthotropic,  cylindrical  shells 
(i . e. ,  the  orthotropic  axes  of  each  layer  coincide  with  the  coordinate 
axes  of  the  shell).  Whitney  and  Sun  [11]  developed  a  shear  deformable 
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theory  for  laminated  cylindrical  shells  that  includes  both  transverse 
shear  deformation  and  transverse  normal  strain  as  well  as  expansional 
strains.  Recently,  Widera  and  Logan  [12,13]  presented  refined  theories 
for  nonhomogeneous  anisotropic  cylindrical  shells. 

As  far  as  the  finite-element  analysis  of  shells  is  concerned, 
layered  composite  shells  have  not  received  nearly  as  much  attention  as 
ordinary  shells.  The  works  of  Dong  [14]  on  statically-loaded  ortho¬ 
tropic  shell  of  revolution,  Dong  and  Selna  [15]  on  free  vibration  of  the 
same,  Wilson  and  Parsons  [16]  on  static  axisymmetric  loading  of  arbi¬ 
trarily  thick  orthotropic  shells  of  revolution,  and  Schmit  and  Monforton 
[17]  on  laminated  anisotropic  cylindrical  shells  are  the  only  ones  that 
considered  the  finite  element  method  before  1 970 ' s  (note  that  the  latter 
reference  is  the  only  one  that  considered  laminated  anisotropic  shells). 

In  1970's  there  was  an  increased  interest  in  the  finite-element  analysis 
of  bending  and  vibration  of  laminated  anisotropic  shells.  Apparently 
the  first  finite-element  application  in  laminated  anisotropic  shells  of 
arbitrary  geometry  is  due  to  Thompson  [18],  who  presented  free  (i.e.,  natural) 
vibration  of  general  laminated  anisotropic  thin  shells.  Other  finite- 
element  analyses  of  layered  anisotropic  composite  shells  include  the 
works  of  Panda  and  Natarajan  [19],  Shivakumar  and  Krishna  Murty  [20], 

Rao  [21],  and  Siede  and  Chang  [28]. 

The  present  paper  is  concerned  with  the  development  and  application 
of  a  shear  deformable  shell  element  for  the  bending  analysis  of  layered 
anisotropic  shells.  The  shear-deformaole  theories  used  herein  are 
generalizations  of  various  classical  shell  theories,  which  are  believed 
to  be  adequate  for  the  prediction  of  overall  response  of  layered 
composite  shells.  In  the  present  investigation,  several  first-order 
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shell  theories  (Sanders,  Love's,  Loo’s,  Morley's  and  Donnell's)  are 
considered  and  transverse  shear  strains  are  accounted  for.  For  certain 
boundary  conditions,  lamination  scheme  and  loading,  exact  form  of  the 
spatial  variation  of  solutions  are  are  obtained  and  are  compared  with 
the  finite-element  solutions. 

GOVERNING  EQUATIONS 

Consider  a  shell  constructed  of  a  finite  number  of  homogeneous, 
uniform  thickness  layers  of  an  orthotropic  material.  Let 
the  x-y  surface  coincide  with  the  lines  of  the  principal  curvature  of 
the  midsurface  of  the  shell,  with  z-axis  normal  (positive  outward)  to 
the  midsurface  of  the  shell  (see  Fig.  1).  It  is  assumed  that  all  of  the 
layers  ir  the  shell  remain  elastic  during  the  deformation,  the  Gener¬ 
alized  Hooke  s  law  is  valid,  and  that  no  slip  occurs  between  any  two 
layers. 

The  displacement  field  in  the  shell  is  assumed  to  be  of  the  form 
u(x,y,z)  =  uQ(x,y)  +  z  *x(x,y) 

v(x.y,z)  =  vQ(x,y)  +  z  fy(x,y)  (1) 

w(x,y,z)  =  wQ(x,y)  +  z  *z(x,y) 

where  u,  v  and  w  are  the  total  displacements  along  x,  y  and  z-direc- 

tions,  respectively;  uQ,  vQ  and  wQ  are  the  displacements  of  a  point  in 

the  coordinate  surface  along  x,  y  and  z  coordinates;  ip  ,  ip  and  ip  are 

x  y  z 

the  bending  rotations.  The  displacement  field  in  eqn.  (1)  can  be 
replaced  by  one  containing  the  higher  order  terms  in  z  (see,  for  example, 
Whitney  and  Sun  [11]).  However,  the  improvement  one  obtains  in  the 
prediction  of  the  transverse  shear  deformation  is  at  the  expense  of 
computational  time  (several  additional  dependent  unknowns  are  introduced 
in  the  higher-order  theories). 
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Strain-displacement  relationships  for  shells  are  very  complex  and  a 
detailed  derivation  cannot  be  attempted  in  a  paper  of  this  type.  References 
to  various  shell  theories  are  made  in  Table  1.  The  strain- 
displacement  relations  for  large-  rotation  (small  strains)  theory  of 
shells  can  be  simplified  to  the  form 

ei  “  ei  +  2<i  »  (i  =  1,2, 4, 5, 6),  (2) 


where 


£1  »T  *  ifj-  *  T  '3 r  • 

=  +  ,3"°.  V 

e2  3y  Rj  2  (3y~  • 


O  3Un  3vr>  3wn  Un  3w„  V^ 

6  3y  3x  3x  +  R^'3y  +  R^’ 


o  .  ,  +  %  '  \ 
4  ■  *y  3y  ■  C1  1^  ’ 


A  ■  *x  *  3T  -  C1  57 


3*y  1  3U« 

r  «  jS  X  ^  X  \  (  Q  t  0  \ 

Ki  jt  +  r^  +  iq-  (jr +  r7}* 

3<’V  1  3V«  W« 

K  =  — Z.  x  —  +  J _ f _ 2.  +  _2.\ 

2  3y  R,  R,  l3y  +  R,'  ’ 


tyy  3<l>  -i  -I  3  V  3urt 

C6  =  37"  +  3^  +  C2  (f£  -  ^  ( JT  '  W~]  ’ 


x4  =  <5  =  0. 


The  underlined  terms  correspond  to  large  rotations.  The  constants  C, 


and  C2  are  shell-theory  tracers  to  be  defined  shortly  (i.e. ,  different 
values  of  C1  trace  different  shell  theories  in  the  literature),  and  R1 
and  R2  are  the  principal  radii  of  curvature  of  the  shell.  The  trans¬ 
verse  normal  strain  e33  is  given  by 

e33  =  *z  +  \  +  (4) 

In  the  present  study  we  make  the  assumption  that  e33  ;  0,  or  equiva¬ 
lently, 

*z  =  '  I^x  +  "'y1 

However,  note  that,  the  transverse  shear  strains  are  accounted  for  in 
the  present  theory.  We  also  assume  that  u^/R1  and  u^/Rg  (for  u^  = 

Vvo>’  wo^Ri  an(^  ^z^Ri  are  sma^  compared  to  the  other  quantities  in  e 
=  1»2,6)  and  are  neglected  in  the  subsequent  equations. 

In  the  absence  of  body  forces  and  moments,  the  equilibrium  equa¬ 
tions  associated  with  a  layered  anisotropic  shell  can  be  simplified  to 
the  form  (consistent  with  the  assumptions  made  above), 

Nl,x  +  N6,y  +  VW  +  C«VC3)/Rr  C2^H6,y  =  °* 

N6,x  +  N2,y  +  VC2/R2}  +  NCrC3)/R2’  C2/R^\x  =  0, 

Ql,x  +  Q2,y  "  Nl/Rl  ‘  N2/R2  *  W(w,Ni,Mi)  =  P, 

Mi,x  +  M6,y-Ql  =0  * 

M6,x  +  M2,y  -  <?2  =  °* 


(5) 


ft 


where  C.  are  the  shell-theory  tracers  given  in  Table  1  and  N.  and  M.  (i 
=  1.2,6)  are  the  stress  and  moment  resultants  defined  in  the  customary 
way  for  a  first-approximation  shell  theory  (z/R..  «  1)  as 
h/2 

(N..M.)  =  /  (1,z)o.  dz  (6) 

-h/2  ' 

where  h  is  the  total  laminate  thickness,  and  ai  are  the  stress  com¬ 
ponents.  The  shear  stress  resultants  are  defined  analogously: 


(Ql » Qo )  =  /  (og.o^Jdz. 

1  c  -h/2  0  4 

The  nonlinear  operator  #(•)  is  defined  by 

.  M,  3w  Mc  3w„ 

.»(«,,») ^ +  c6  -  4^ 


M0 

0 


M r  3W 
0\  0- 


♦  ly  f<N2  -  4>5T  *  (N6  - 


Table  1  List  of  shell-theory  tracers 


Theory 


Sanders  [24]  1.0  0.5  1.0  1.5 

Love’s  1st  Approximation  [25]  1.0  0.0  1.0  1.0 

Loo's  Approximation  [26]  0.0  0.0  1.0  1.0 

Mor1 ey  [27]  0.0  0.0  1.0  0.0 

Donnell  [28]  0.0  0.0  0.0  0.0 

The  thermoelastic  constitutive  equations  for  a  layered  anisotropic 
shell  can  be  derived  using  the  constitutive  equations  of  individual 
layers  (in  much  the  same  way  as  for  plates;  see,  Reddy  [8]).  The 
equations  are  given  by 
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jQU  .  fA55  S45]  (c4) 

*  Q2*  L  A45  A44  J  (  e5J 

iy- 


where  the  stiffness  coefficients  A..,  B. .  and  D. .  are  given  by 

•  J  *  J  *  J 

h/2  7 

<V  V  °ij)  *  /.h/2  OljO'4.2  )<!*..  ('.j  •  '.2’6>. 


Ai j  ■  KZ  /  Q-n  dz  ,  (i,j  =  4,5), 

1J  -h/2  1J 

and  the  thermoelastic  stress  resultants  and  stress  couples  are  defined 


(Nj,Mj)  «  /  (l,z)Q..a.  dz.  (11) 

l  1  -h/2  10  0 

Here  a.  denotes  the  thermal  expansion  coefficients  and  Q. .  denote  the 
J  .  v 

plane-stress-reduced  stiffnesses  for  individual  layers.  Since  Q^j  are 
different  for  different  layers,  expressions  in  eqns.  (9)  and  (10)  should 
be  integrated  in  parts,  for  example, 

„  m  zt+l  r„\ 


Aij  =  /  Q  j  dz 

1J  i-l  z.  1J 


where  m  is  the  total  number  of  layers  in  the  shell,  and  K  is  the  shear 
correction  coefficient. 

Equations  (2),  (5),  and  (9)  completely  describe  the  equilibrium 
equations  of  a  layered  anisotropic  shell.  These  equations  cannot  be 
solved  in  an  exact  form  except  for  the  small-displacement  theory,  simple- 
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supported  boundary  conditions,  and  sinusoidal  distribution  of  the  trans¬ 
verse  load  and  thermal  loads  (see  Hsu,  Reddy  and  Bert  [23]).  In  order 
to  solve  practically  important  problems  that  involve  complex  loadings 
and  geometry,  one  must  consider  approximate  methods  of  analysis.  In  the 
next  section,  a  finite-element  formulation  is  presented  for  the  equations 
described  in  this  section. 

FINITE-ELEMENT  FORMULATION 

In  the  interest  of  brevity,  we  omit  the  algebraic  details  and 
describe  briefly  the  finite-element  model  developed  herein.  The  derivation 
proceeds  along  the  same  lines  as  described  for  layered  anisotropic 
plates  in  [29].  Consider  a  finite  element  analog  of  the  midsurface  R  of 
the  shell.  In  a  typical  (finite)  element  Re  of  the  mesh,  the  generalized 
displacements  u0»v0»w0i  ^  and  ^  are  interpolated  by  expressions  of  the 
form, 


uo  =  l.  uoi*i  »  vo  '  J  voiV  etc'  03) 

where  uQ^  is  the  value  of  uQ  at  node  i  of  the  element  Re,  <f>.  is  the 
■finite-element  interpolation  function  associated  with  node  i,  and  n  is 
the  total  number  of  nodes  in  the  element.  Here  we  have  employed  the 
same  interpolation  functions  for  each  of  the  dependent  variables. 
Substituting  the  expressions  (13)  ‘nto  the  governing  equations  (5), 
multiplying  each  equation  by  <*>.,  and  setting  the  integral  of  the  result 
over  Re  to  zero  (Galerkin  integrals),  we  obtain  (after  carrying  integration 
by  parts  of  certain  terms  to  relax  the  continuity  on  $.)  the  following 
form  of  element  equations: 


i 


04) 


CK]{A>  *  {F}. 

Here  denotes  the  stiffness  matrix  (5n  by  5n),  {A}  «  *u0i ,voi ,woi ’ 
Vi’^yv  '  an<*  the  vector  °f  generalized  forces. 

In  the  present  study  the  four-node,  eight-node  and  nine-node 
rectangular  isoparametric  elements  were  employed.  Analogous  to  the 
shear  deformable  theory  of  layered  composite  plates,  the  present  theory 
can  be  viewed  as  a  shear  deformable  theory  derived  from  the  classical 
shell  theory  by  treating  the  slope-deflection  relations  (e^  and  eg  =  0) 
as  constraints,  and  including  the  constraints  into  the  variational 
formulation  of  the  shell  equations  by  the  penalty  function  method  (see 
Reddy  [29]).  The  elements  derived  using  such  theory  are  very  stiff  for 
thin  shells,  but  yield  good  results  for  moderately  thick  shells.  To 
overcome  the  difficulty,  the  so-called  reduced  integration  technique 
must  be  employed  in  the  evaluation  of  the  stiffness  coefficients  asso¬ 
ciated  with  the  shear  energy  terms  (i.e. ,  penalty  terms).  More  speci¬ 
fically,  the  lxl  Gauss  rule  must  be  used  for  shear  energy  terms  and  the 
standard  2x2  Gauss  rule  must  be  used  for  the  bending  terms  when  the 
four-node  linear  isoparametric  element  is  used. 

NUMERICAL  RESULTS 

All  of  the  numerical  results  presented  herein  were  obtained  on  an 
IBM  370/158  computer  using  the  double  precision  arithmatics.  Whenever 
biaxial  symmetry  existed  in  the  problem,  only  one  quadrant  of  the  shell 
was  modeled.  Numerical  results  are  presented  using  the  linear  and 
nonlinear  theory.  In  the  nonlinear  analysis,  a  direct  iteration  method 
was  used.  The  iteration  begins  with  zero  solution  vector  (so  that  at 
the  end  of  the  first  iteration  the  linear  solution  is  obtained)  and 
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computes  the  linear  element  matrices.  The  solution  at  the  end  of  the 
first  iteration  is  used  to  recompute  the  element  stiffness  matrices  (for 
the  same  load  step)  for  the  next  iteration.  This  procedure  Is  continued 
until  the  difference  between  two  consecutive  solution  vectors,  for  a 
given  load,  differ  by  one  percent.  A  shear  correction  factor  of  K  = 

5/6  is  used  In  all  of  the  cases,  and  only  one  quadrant  of  the  shell  is 
modeled  (whenever  the  biaxial  symmetry  existed). 

To  assess  the  numerical  accuracy  of  the  present  model  with  respect 
to  the  element  type  (i.e.,  linear  and  quadratic  elements),  mesh  (L2  =  2 
by  2  mesh  of  the  linear  elements,  Q2  =  2  by  2  mesh  of  the  eight-node 
quadratic  elements,  and  9  =  nine-node  elements),  and  integration  (F  = 
full  integration,  R  =  reduced  integration  on  all  terms),  several  numerical 
experiments  were  carried  using  the  Sanders  shell-theory,  and  the  results 


are  listed  in  Table  2.  The  following  nondimensionalizations  are  used: 


_  wE,h  , 
w  =  — V  x  lcr 


Table  2  contains  the  nondimensionalized  linear  deflection  and  stresses 


for  a  four-layer  (0°/90°/  90°/0‘'),  equal  thickness,  spherical  shell  (R^ 
=  Rg  =  R,  a/h  3  10)  of  material  1  (G12  =  v^2  =  v13): 

material  1:  E-j/E2  8  25,  G12/E2  -  0.5, 

G2j/E2  =  0.2,  v-j2  ~  0.25.  (15) 

The  shell  is  subjected  to  sinusoidal  distribution  of  normal  loading, 


P  =  p.  sin  ~  si 

U  a 


«»?■ 


where  a  and  b  are  the  dimensions  shown  in  Fig.  1.  The  following  boundary 


conditions  were  imposed: 
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SS-1 : 


natural : 


essential 


=  0  at  x  =  0,a 

Ng  =  M2  =  0  at  y  =  0,b 

wo  =  vo  3  0  at  x  =  0,a 

wQ  =  uQ  =  0  at  y  »  0,b 


(18) 


For  this  special  case  of  boundary  conditions,  loading  and  cross-ply 
construction,  one  can  obtain  the  exact  form  of  the  solution;  see  Hsu, 
Reddy  and  Bert  [23]. 


From  Table  2  one  can  conclude  that  the  reduced  integration  improves 
the  deflections  as  well  as  stresses  for  thin  shells,  when  coarse  meshes 
and  lower  order  (i.e. ,  linear)  elements  are  used.  For  quadratic  elements 
the  reduced  integration  has  only  little  effect  even  for  a/h  =  100. 
Numerical  convergence  of  the  element  is  apparent  from  an  inspection  of 
the  deflections  and  stresses  in  Table  2.  The  agreement  between  the 
present  finite  element  results  and  the  exact  solution  is  gratifyingly 
close.  It  should  be  pointed  out  that  both  the  exact  solution  and  the 
finite-element  solutions  are  obtained  for  the  same  shell  theories.  The 


subsequent  results  were  obtained  using  the  reduced  integration  technique. 

Next,  the  effect  of  thickness  and  type  of  shell  theory  on  linear 
deflections  and  stresses  of  a  four-layer  (0°/90  790°/0°,  material  1) 
cylindrical  shell  (R-j  =  R,  Rg  =  1030)  subjected  to  sinusoidal  distribution 
of  normal  force  and  simply  supported  (SS-1)  boundary  conditions  is 
investigated  (see  Table  3).  From  the  results  (obtained  using  2x2Q)  of 
Table  3  it  is  clear  that  for  shallow  shells  (R/h  =  100)  the  difference 
between  various  theories  is  negligibly  small  for  the  radius-to-thickness 
ratio  considered.  For  radius-to*thickness  ratio  of  10,  the  difference 
between  the  deflections  obtained  by  various  shell  theories  is  noticeable. 
In  the  examples  to  follow,  the  Sanders  shell  theory  was  used. 
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Table  3  Dimensionless  Deflections  and  Stresses  for  4-Layer  Cross-Ply  (0°/90o/90o/0°) 
Freely  Supported  Cylindrical  Shells  Under  Sinusoidal  Loading  by  Different 
Shell  Theories  (material  1,  a/b  =  1,  a/h  =  10). 
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Table  4  shows  the  effect  of  number  of  layers,  orientation  of 
layers,  and  radius-to-thickness  rat’o  on  the  nondimensionalized  deflection 
of  a  simply  supported  (SS-1),  angle-ply  (8/-Q/+/-- . . )  spherical  shell 
(R.|  *  R2  =  R)  under  sinusoidal  loading.  The  shell  is  constructed  of 
material  2  (G^  =  G13‘  v12=v13): 

material  2:  E1/E2  =  40,  G12/E2  =  0.6,  . 

G23/e2=  0.5,  V12  =  0.25. 

The  finite-element  solution  differ  noticeably  from  the  exact  solution 
(Fig.  2)  for  small  number  of  layers  (the  difference  decreases  with 
increasing  number  of  layers)  and  large  values  of  the  angle,  e.  This 
error  is  due  to  the  fact  that  the  exact  solution  given  in  [23] 
satisfies  different  symmetry  conditions  than  those  imposed  in  the 
finite  element  method.  This  difference  is  directly  proportional  to 
stiffness  coefficient  B^,  whose  magnitude  is  inversely  proportional 
to  the  number  of  layers  and  directly  proportional  to  e(0  £  e  <_ 45°). 

Figure  3  shows  the  effect  of  aspect  ratio  (b/a)  and  radius-to- 
thickuess  ratio  (R/h)  on  the  nondimensional  deflections  (wt  =  wh/a-jT^) 
for  simply  supoorted  (SS-1)  spherical  shell  (R1  =  R2  =  R,  material  1) 
under  simusoidal-thermal  loading. 

T  =  T0  sin  §*■  sin  (20) 

It  is  clear  from  the  figure  that  thickness  shear  deformation  has  a 
pronounced  effect  on  the  nondimensionalized  deflection,  w^;  the  effect 
decreases  with  increasing  number  of  layers  (the  two-layer  solution 
forms  the  upper  bound  while  the  single-layer  solution  forms  the  lower 
bound). 


Table  4.  Nondimensional  deflection  as  a  function  of  number  of  layers, 
angles,  and  radius-to- thickness  ratio  for  freely  supported 
angle-ply  ( 9/- e/. ..)  spherical  shell  (a/b  =  1,  a/h  =  10, 
material  2). 


lC2E2h3 


R/h  Source 


0  =  5' 
N=4 


6  =  30' 
n=4 


0  =  45' 
n=4 


5 

CFS 

0.1182 

0.1182 

0.1182 

0.03165 

0.03165 

0.03165 

0.01183 

0.01183 

0.01183 

FEM 

0.1424 

0.1255 

0.1194 

0.03509 

0.03147 

0.03140 

0.009382  0.01079 

0.01138 

10 

CFS 

0.2673 

0.2673 

0.2673 

0.09552 

0.09552 

0.09552 

0.04170  0.04170 

0.04170 

FEM 

0.3436 

0.2924 

0.2720 

0.1306 

0.1002 

0.09566 

0.05215 

0.04061 

0.04079 

20 

CFS 

0.3820 

0.3820 

0.3820 

0.1908 

0.1908 

0.1908 

0.1127 

0.1127 

0.1127 

FEM 

0.4715 

0.4112 

0.3873 

0.3089 

0.2070 

0.1920 

0.1783 

0.1175 

0.1118 

30 

CFS 

0.4146 

0.4146 

0.4146 

0.2339 

0.2339 

0.2339 

0.1646 

0.1646 

0.1646 

FEM 

0.4898 

0.4400 

0.4189 

0.4038 

0.2562 

0.2354 

0.2825 

0.1753 

0.1640 

40 

CFS 

0.4273 

0.4273 

0.2540 

0.2540 

0.2540 

0.2540 

0.1962 

0.1962 

0.1962 

FEM 

0.4953 

0.^496 

0.4308 

0.4505 

0.2791 

0.2556 

0.3507 

0.2110 

0.1960 

50 

CFS 

0.4334 

0.4334 

0.4334 

0.2645 

0.2645 

0.2645 

0.2153 

0.2153 

0.2153 

FEM 

0.4989 

0.4535 

0.4365 

0.4752 

0.2910 

0.2661 

0.3940 

0.2329 

0.2154 

60 

CFS 

0.4369 

0.4369 

0.4369 

0.2706 

0.2706 

0.2706 

0.2274 

0.2274 

0.2274 

FEM 

0.5003 

0.4553 

0.4395 

0.4895 

0.2979 

0.2722 

0.4223 

0.2468 

0.2276 

60 

CFS 

0.4403 

0.4403 

0.4403 

0.2769 

0.2769 

0.2769 

0.2408 

0.2408 

0.2408 

FEM 

0.5008 

0.4567 

0.4424 

0.5039 

0.3049 

0.2786 

0.45F0 

0.2626 

0.2413 

100 

CFS 

0.4419 

0.4419 

0.4419 

0.2799 

0.2799 

0.2799 

0.2476 

0.2476 

0.2476 

FEM 

0.5016 

0.4570 

0.4437 

0.5106 

0.3082 

0.2815 

0.4722 

0.2707 

0.2483 

00 

CFS 

0.4448 

0.4448 

0.4448 

0.2355 

0.2855 

0.2855 

0.2505 

0.2605 

0.2605 

(plate) 

FEM 

0.5019 

0.4546 

0.4454 

0.5199 

0.3137 

0.2870 

0.5119 

0.2879 

0.2620 

^-FEM  (R/h=100,  N=2) 

- - - 

- - - 

P  FEM  (R/h=100,  N=16) 

FS  (R/h=100,N=2) 


Oj  _°FEM(R./h=10,  N=16) 

.  CFS(R/h=10,  N=2)  - 


FEM  (R/h=10,  N=2) 


e°-  5  10  15  20  25  30  35 


40  45  50 


Figure  2.  Effect  of  lamination  scheme  and  radius-to-thickness  ratio 
on  the  nondimensionalized  deflection  of  a  sperical  shell 

SSriST  transverse  fading  (SS-1,  aP/h  '  io! 


(b)  Deflection  versus  aspect  ratio. 

Figure  3.  Effect  of  thickness  and  aspect  ratio  on  the  deflection 
of  cross-ply  doubly-curved  shells  under  thermal  loading 
(material  1,  R, .  »  R„  *  R) 
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The  natural  vibration  of  layered  composite  shells  was  also  investi¬ 
gated  using  the  element  presented  herein.  Tables  5  and  6  contain  the 
values  of  the  nondimensional ized  fundamental  frequencies  of  doubly  curved 
(spherical)  shells  of  cross-ply  and  angle-ply  construction,  respectively. 
Note  that  as  the  number  of  layers  is  increased,  the  frequencies  do  not 
change  appreciably  (for  the  same  total  thickness  o'  che  shell).  The 
nondimensional ized  frequencies  of  the  cross-ply  shells  are  about  half  of 
those  of  the  angle-ply  shells. 


Table  5. 


Effect  of  radius-to-thickness 
fundamental  frequency,  \  =  „r 
shells  t(Rj=R2=R,  material  2, 


ratio  on  the  dimensionless 
p/E?  ,  of  cross-ply  spherical 
a/b  =  1,  a/h  =  10). 


R/  h  Source  0°/90°  0°/90o/0°  0°/90#/90°/0# 


5 

CFS 

1.0536 

1.0739 

1 .0757 

FEM 

1.0538 

1.0746 

1.0713 

10 

CFS 

1.4452 

1.6085 

1 .61*6 

FEM 

1.4442 

1.6128 

1.6167 

20 

CFS 

2.1484 

2.675 

2.688 

FEM 

2.1518 

2.6818 

2.6812 

30 

CFS 

2.934 

3.8184 

3.8379 

FEM 

2.9352 

3.8175 

3.8073 

40 

CFS 

3.7636 

4.9944 

5.0204 

FEM 

3.7538 

4.9824 

5.0220 

50 

CFS 

4.6151 

6.186 

6.218 

FEM 

4.6035 

6.0045 

6.2695 

60 

CFS 

5. 1786 

7.3854 

7.4238 

FEM 

5.4592 

7.3926 

7.5222 

80 

CFS 

7.2247 

9.7960 

9.9480 

FEM 

7.1954 

9.8280 

9.8472 

100 

CFS 

8.9461 

1 2.225 

12.274 

FEM 

8.9001 

12.163 

12.238 

freely-supported  boundary  conditions  are  used. 


Table  6.  Effect  of  radius-to-thickness  ratio  on  the  nondimen- 
sional  fundamental  frequencies  of  angle-ply  spherical 
shel  1  s+(R-|  =R2=R»  material  2,  a/b  =  1 ,  a/h  =  10). 


R/h 

Source 

45*7-45° 

45°/-450/45°/-450 

45*7-45°/ 

n=8 

5 

CFS 

2.2215 

FEM 

2.2229 

2.2232 

2.2236 

10 

CFS 

4.4429 

FEM 

3.6286 

4.4456 

4.4461 

20 

CFS 

6.0540 

FEM 

4.5238 

5.7334 

5.9716 

30 

CFS 

7.6344 

FEM 

5.5908 

7.2330 

7.536 

40 

CFS 

9.4032 

FEM 

6.81 00 

8.9036 

9.2804 

50 

CFS 

11.275 

FEM 

8.100 

10,6895 

11.1315 

60 

CFS 

13.266 

FEM 

9.4794 

12.516 

13.0434 

80 

CFS 

17.1688 

FEM 

12.2704 

16.2752 

16.9608 

100 

CFS 

21.202- 

FEM 

15.129 

20.095 

20.948 

^  Freely-supported  boundary  conditions  are  used. 
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SUMMARY  AND  CONCLUSIONS 

A  finite  element  based  on  a  shear  deformation  theory  of  layered, 
anisotropic,  composite  shells  is  presented.  The  element  contains  as 
special  cases  Generalized  first-approximation  shell  theories  of  love, 
Sanders,  Donnell,  Mori ey,  and  Loo.  The  generalization  is  to  include 
the  transverse  shear  strains.  Numerical  results  are  presented  for 
cross-ply  as  well  as  angle-ply  shells  under  various  mechanical  and 
thermal  loadings.  It  Is  concluded  that  the  effect  of  shear  deforma¬ 
tion  on  the  deflection  and  stresses  is  significant,  and  that  reduced 
integration  improves  the  deflections  and  stresses  for  thin  shells  when 
coarse  meshes  and  lower  order  elements  are  employed.  Use  of  higher- 
order  theories  and  extension  of  the  present  element  to  the  nonlinear 
(geometrically)  analysis  of  shells  is  far  from  complete.  Although  the 
nonlinear  terms  in  the  sense  of  von  Karman  are  included  in  the  formulation 
presented  herein,  no  numerical  results  are  included.  These  will  appear 
elsewhere  when  the  investigation  is  completed. 
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